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aluminum  alloy  2024-0,  under  a central  concentrated  load,  with  large  de- 
flection. The  load  was  provided  by  a small  diameter  hard  steel  rod.  Mea- 
surements were  made  of  load,  deflections,  and  strains,  and  membrane  and 
bending  strains  were  calculated  from  the  test  data.  Twenty  electrically 
bonded  strain  gauges  were  used  to  measure  the  distribution  of  radial  and 
circumferential  strain  components  along  a radial  line  on  both  faces  of 
each  plate.  The  test  data  are  presented  in  comparison  with  theoretical 
predictions  generated  by  the  Grumman-developed  finite  element  computer 
code  PLANS.  The  correlation  between  theory  and  test  data  was  excellent 
for  deflections,  and  generally  good  for  strains,  using  the  circular  line 
load  predictions.  The  uniform  pressure  case  slightly  overpredicted  the 
central  deflections  and  greatly  overpredicted  the  peak  strains.  The 
greatest  deviation  between  predicted  and  measured  strains  was  in  the  re- 
gion of  the  concentrated  load,  where  local  transverse  normal  and  shear 
stresses  are  not  included  in  the  theory,  and  where  the  errors  caused  by 
the  assumed  load  distribution  would  be  maximum.  The  plates  exhibited 
initial  loss  of  stiffness  under  the  plastic  bending  behavior,  followed  t 
a rapidly  increasing  membrane  action  resulting  from  large  deflections, 
which  provided  much  additional  resistance  to  the  applied  load. 
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Tests  were  performed  on  three  simply  supported  circular  plates 
of  aluminum  alloy  2024-0,  under  a central  concentrated  load,  with 
large  deflection.  The  load  was  provided  by  a small  diameter  hard 
steel  rod.  Two  plates,  having  a diameter- to-thickness  ratio  of 
D/h  = 40.6,  were  loaded  to  a central  deflection  of  2.6  times  the 
thickness,  and  one  plate  of  D/h  = 20.0  was  loaded  to  a central 
deflection  of  1.6  times  the  thickness.  Measurements  were  made  of 
load,  deflections,  and  strains,  and  membrane  and  bending  strains 
were  calculated  from  the  test  data.  Twenty  electrically  bonded 
strain  gauges  were  used  to  measure  the  distribution  of  radial  and 
circumferential  strain  components  along  a radial  line  on  both  faces 
of  each  plate.  The  test  data  are  presented  in  comparison  with 
theoretical  predictions  generated  by  the  Grumman -developed  finite 
element  computer  code  PLANS.  The  theoretical  model  used  two  dif- 
ferent assumed  distributions  of  contact  pressure  under  the  loading 
rod  to  approximate  the  test  load  distribution.  They  were  a uniform 
contact  pressure  and  a concentrated  line  load  along  the  rod  edge 
circle.  The  correlation  between  theory  and  test  data  was  excellent 
for  deflections,  and  generally  good  for  strains,  using  the  circular 
line  load  predictions.  The  uniform  pressure  case  slightly  overpre- 
dicted the  central  deflections  and  greatly  overpredicted  the  peak 
strains.  The  greatest  deviation  between  predicted  and  measured 
strains  was  in  the  region  of  the  concentrated  load,  where  local 
transverse  normal  and  shear  stresses  are  not  included  in  the  theory, 
and  where  the  errors  caused  by  the  assumed  load  distribution  would 
be  maximum.  The  plates  exhibited  initial  loss  of  stiffness  under 
the  plastic  bending  behavior,  followed  by  a rapidly  increasing  mem- 
brane action  resulting  from  large  deflections,  which  provided  much 
additional  resistance  to  the  applied  load. 
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1.  INTRODUCTION 


A series  of  tests  were  performed  for  the  purpose  of  exhibit- 
ing details  of  the  plastic  bending  characteristics  of  a representa- 
tive structural  element  and  to  evaluate  the  predictive  capability 
of  the  AXSHEL  (Axi symmetric  Shell  Analysis)  module  of  the  Grumman- 
developed  PLANS  (Plastic  Large  Deflection  Analysis  of  Structures) 
system  of  computer  codes. 

The  simply  supported  circular  plate  under  a central  concen- 
trated load  was  chosen  for  these  purposes  because  it  is  one  of  the 
simplest  structural  elements  exhibiting  a hiaxial  stress  state  in 
bending.  The  simple  loading  and  geometry  would  allow  for  easy 
fabrication  and  fixturing,  while  the  load  concentration  would  pro- 
vide a severe  test  of  the  theory  because  of  the  large  strain  gra- 
dients and  contact  stresses.  The  material  chosen  was  2024-0 
(annealed)  aluminum  alloy,  which  has  a very  small  elastic  strain 
range  and  a large  plastic  strain  before  failure,  thereby  allowing 
large  plastic  strains  to  be  developed  at  relatively  low  loads.  It 
has  a gradually  changing  tangent  modulus  in  the  plastic  strain 
range,  which  is  more  typical  for  aerospace-type  materials  than  the 
abrupt  yielding  with  zero  plastic  tangent  modulus  (to  2 or  3 per- 
cent strain)  exhibited  by  the  mild  steel  used  in  most  previously 
reported  tests. 

The  tests  were  required  to  provide  data  on  the  distribution 
of  strain  components  and  the  deflections.  This  is  because  the  de- 
flections reflect  an  integrated  behavior  of  all  points  in  the 
structure, and  the  transverse  motion  of  any  local  point  is  there- 
fore a somewhat  gross  measure  of  the  structural  response.  The  local 
strains,  however,  are  much  more  dependent  on  purely  local  behavior, 
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and  thus  examination  of  the  strain  distribution  over  the  structure 
can  provide  more  information  on  the  variation  of  structural  be- 
havior from  point  to  point. 

While  there  has  been  a moderate  amount  of  test  data  available 
to  verify  various  elastic-plastic  analysis  methods  (Refs.  1 through 
19),  most  reports  show  only  deflections  but  do  not  report  strain 
distributions.  The  notable  exceptions  are  tests  by  Ohashi  and 
Murakami  (Refs.  13  and  14),  Ohashi  and  Kawashima  (Ref.  16),  and 
May  (Ref.  15)  for  the  moderate  thickness  range  of  mild  steel  plates. 
Cooper  and  Shifrin  (Ref.  6)  reported  strains  and  curvatures  only 
near  the  edge  of  very  thin  mild  steel  plates,  for  which  the  mem- 
brane behavior  was  dominant.  Some  earlier  documents  (Refs.  2 and 
3)  by  Ramberg,  McPherson,  and  Levy  reported  residual  deformation 
and  strain  data  on  very  thin  aluminum,  magnesium,  and  stainless 
steel  plates.  They  applied  the  load  cyclically,  increasing  the  peak 
load  after  each  release,  so  that  the  usual  monotonically  increasing 
stress-strain  data  does  not  apply.  The  flow  patterns  are  of  addi- 
tional interest, as  reported  by  Lance  and  Onat  (Ref. 12).  Recently, 
Durelli,  Parks,  and  Chen  (Ref.  19)  have  reported  deformations  and 
strains  for  a centrally  loaded  thick  plate  (D/h  = 6).  They  used 
the  Moire  optical  strain  measuring  method,  as  well  as  the  usual 
bonded  resistance  strain  gauges.  However,  their  large  plate  thick- 
ness would  make  the  shear  stresses  too  prominent  for  comparison 
with  the  AXSHEL  bending  and  membrane  theory.  One  of  the  most  in- 
teresting tests  on  mild  steel  was  that  of  Sherbourne  and  Srivastava 
(Ref.  18)  who  carried  the  load  level  high  enough  to  exhibit  initial 
plastic  softening  under  bending  action,  restiffening  under  increas- 
ing membrane  action,  and  a second  plastic  softening.  But,  unfor- 
tunately, they  did  not  report  any  strain  data. 
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In  summary,  no  test  data  were  found  for  both  deflections  and 
strains  on  plates  of  a gradually  yielding,  work  hardening  material 
exhibiting  combined  bending  and  membrane  action.  The  available 
teported  test  data  were  considered  inadequate  for  the  stated  pur- 
poses, so  a related  test  program  was  undertaken,  as  follows. 
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2 . PROCEDURE 


Test  Specimens 

Experiments  were  performed  on  three  5.35-inch  diameter  flat 
circulax  plate  specimens,  two  of  which  were  ^-inch  nominal  thick- 
ness and  one  of  ^--inch  nominal  thickness.  The  material  stress- 
strain  curves  were  determined  by  tests  on  16  tension  coupons,  com- 
posed of  8 coupons  for  each  thickness,  and  two  compression  coupons 
for  the  |-inch  plates. 

It  was  decided  to  make  the  specimens  from  aluminum  alloy  2024-0 
(annealed),  because  its  low  yield  strength  at  room  temperature  would 
produce  large  plastic  strains  at  small  loads.  Also  important  was 
its  reported  isotropy  and  nearly  equal  tension  and  compression 
properties.  After  examining  a number  of  candidate  4 by  8 foot- 
plates, plates  of  f-inch  and  of  £-inch  nominal  thickness  were 
selected  with  clean,  unmarred  surfaces,  and  a 2-foot  square  section 
was  cut  from  each.  Four  circular  plate  blanks,  and  12  tension  cou- 
pon blanks  were  saw-cut  from  this  square  plate,  following  the  pat- 
tern shown  in  Fig.  1. 

The  plate  specimens  and  tension  coupons  were  machined  on  their 
edges  only,  leaving  the  original  face  surfaces  intact.  The  plate 
specimens  were  finished  to  a nominal  5.35-inch  diameter.  The 
final  tension  coupon  dimensions  are  shown  in  Fig.  2.  The  plates 
and  tension  coupons  were  marked  on  the  edges  to  keep  track  of  cheir 
numbers  according  to  the  layout  of  Fig.  1.  The  scrap  material  re- 
maining from  the  original  plates  was  saved,  and  was  used  later  in 
the  project  to  make  and  ^-inch  diameter  by  1-inch  long  cylin- 

drical compression  coupons. 
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In  summary,  no  test  data  were  found  for  both  deflections  and 
strains  on  plates  of  a gradually  yielding,  work  hardening  material 
exhibiting  combined  bending  and  membrane  action.  The  available 
reported  test  data  were  considered  inadequate  for  the  stated  pur- 
poses, so  a related  test  program  was  undertaken,  as  follows. 
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Although  a total  of  8 plate  blanks  were  cut  from  the 
and  £-inch  plate  stock,  only  four  circular  plate  specimens, 
numbered  1 and  4 of  each  thickness,  were  finish  machined  and 
instrumented  at  this  time.  Of  these  four,  specimen  number  1A250 
(number  1 of  0.250-inch  nominal  thickness)  was  not  tested  at 
this  time.  Thus,  only  three  plate  specimens,  composed  of  two 
5-inch  plates  (1A125  and  4A125)  and  one  £-inch  plate  (4a250), 
are  reported  on.  Only  those  tension  coupon  blanks  immediately 
adjacent  to  plate  blanks  1 and  4 were  machined  into  finished 
tension  coupons.  These  are  the  8 blanks  of  each  thickness 
numbered  1,  2,  3,  6,  8,  9,  11,  and  12  in  Fig.  1. 

Plate  Tests 

After  each  circular  plate  specimen  was  finished,  it  was  care- 
fully measured  to  determine  its  diameter,  thickness  variation,  and 
deviation  from  flat.  These  measurement  data  are  listed  in  Table  1, 
for  the  three  plates  reported  here. 

Tne  plate  test  arrangement  is  shown  in  Figs.  3,  4,  and  5.  The 
plate  specimen  (1),  rested  on  the  hard  steel  support  fixture  (2), 
which  was  installed  on  the  lower  crossbeam  (3)  of  the  testing 
machine  frame.  The  testing  machine  was  an  Instron  Model  TTDL,  a 
screw-driven,  constant  crosshead  speed  machine  with  a 20,000  pound 
load  capacity.  The  loading  rod  (4)  was  attached  to  the  load  cell 
(5),  which  was  mounted  to  the  underside  of  the  moving  crossbeam  (t) , 
(see  Fig.  4).  The  testing  machine  control  console  (7)  also  con- 
tained the  load  cell  conditioner  and  recording  chart.  The  cabinet 
on  the  right  housed  the  strain  gauge  conditioning  and  recording 
equipment,  including  a 20-channel  scanner/indicator  (8),  and  a 
digital  printer  (9).  To  the  left  of  the  plate  specimen  were  the 
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Linear  Variable  Differential  Transformer  (LVDT)  deflection  trans- 
ducer conditioning  and  recording  units,  including  the  dc  power 
supplies  (10) , a signal  balancing  unit  (11) , and  a strip-chart 
recorder  (12) . 

A closeup  view  appears  in  Fig.  5 of  the  plate  support  fixture, 
the  three  LVDT  displacement  transducers,  and  the  loading  rod  at- 
tached to  the  load  cell.  These  are  shown  in  position  prior  to  the 
installation  of  the  plate  specimen.  The  fixture  opening  was 
5.22-inches  in  diameter,  supporting  cne  plate  specimens  of 
5.35-inch  nominal  diameter.  The  harlened  steel  loading  rod  was 
initially  0.430-inch  in  diameter  for  plate  No.  1A125,  and  was 
then  reduced  to  0.375  inch  for  the  remaining  two  plates. 

The  instrumentation  used  in  the  plate  tests  were  3 LVDT  de- 
flection transducers,  21  surface -bonded  strain  gauges,  and  one 
load  cell. 

The  three  LVDT's  are  shown  in  place  in  Fig.  5,  and  were  the 
dc  in-dc  out  type  (Tresco  Inc.)  with  built-in  oscillator  and 
rectifier,  a 0.10-inch  linear  range,  and  a 0.20-inch  usable 
range.  The  LVDT's  were  located  on  a radial  line  at  nominal  posi- 
tions of  0,1,  and  2 inches  from  the  center  of  the  support 
fixture.  The  contact  tip  on  the  probe  of  the  centerline  LVDT  was 
flat,  to  pick  up  the  peak  deflection  of  the  center  of  the  plate, 
while  the  other  two  probe  tips  were  sharp-pointed,  to  remain  in 
constant  touch  with  their  initial  contact  points  on  the  plate 
bottom.  The  LVDT's  were  supplied  with  24  volts  regulated  dc 
power;  the  output  signals  were  sent  to  a potentiometer  box  that 
balanced  the  starting  signals  to  zero  and  then  sent  to  the  strip 
chart  recorder. 
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Motion  of  the  load  rod  tip  was  also  used  as  a source  of  de- 
flection data.  This  was  done  by  relating  the  known  constant 
crossbeam  speed  to  the  known  constant  speed  of  the  Instron  chart 
recording  the  load,  and  then  adjusting  for  the  calibrated  deflec- 
tion versus  load  characteristics  for  the  machine  in  that  particu- 
lar configuration. 

The  strain  gauge  layout  for  the  plates  is  shown  in  Fig.  6 and 
a typical'  gauge  installation  is  shown  in  Figs.  7 and  8.  The  upper 
surface  was  defined  as  the  surface  contacted  by  the  loading  rod. 
Gauges  2 through  17  were  90  degree  dual-element  type  (FAET-12D- 
12S13ET,  BLH  Elect! onics  Inc.)  and  gauges  1 and  18-21  were  single- 
element type  (FAE-12S-12S13ET,  BLH  Electronics  Inc.),  all  of  the 
same  lot  number  (A270).  These  are  foil  gauges  sandwiched  between 
polyimide  carrier  and  cover  strips,  with  factory-attached  copper 
terminals,  a f-inch  square  sensing  grid,  a 120-ohm  nominal  re- 
sistance, and  a thermal  expansion  rate  of  13.0  x 10~  per  °F. 

The  manufacturer's  data  show  a thermally  induced  strain  error  of 
-10  microstrain  or  less  on  aluminum,  and  a gauge  factor  change  of 
±0.2  percent  or  less,  both  between  50  and  125°F. 

This  type  of  gauge  configuration  was  selected  after  an  ex- 
ploratory test  series,  described  in  Appendix  A.  These  tests  showed 
that  such  gauges  could  survive  almost  as  much  peak  strain  as  the 
traditional  handmade  wire  jumper  gauge  configuration,  with  much 
less  installation  effort. 

The  bridge  circuit  was  tho  single  active  arm  type  with  two 
lead  wires,  with  the  nonlinearity  effect  compensated  by  data  pro- 
cessing, The  bridge  excitation  was  5 volts.  The  strain  gauges 
were  placed  so  that  both  the  radially  and  circumferentially  oriented 
gauges  of  a pair  would  be  located  at  the  same  radius.  This  caused 
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the  gauges  to  be  rotated  slightly  with  respect  to  the  radial  and 
circumferential  axes  through  the  gauge  centers.  The  effect  of 
this  was  corrected  in  the  data  processing.  The  dual-gauge  sets 
provided  for  strain  measurements  mutually  perpendicular  to  each 
gauge,  thus  permitting  transverse  gauge  error  to  be  corrected  in 
data  processing.  The  various  gauge  correction  procedures  are  de- 
scribed in  Appendix  B. 

Specimen  surfaces  were  roughened  lightly  by  an  abrasive  powder 
blast  in  the  areas  to  receive  the  strain  gauges,  with  an  associated 
thickness  reduction  of  less  than  0.0002  inch.  The  gauges  were 
then  bonded  to  the  surface  using  a two-part  epoxy  cement  (GA-2, 
Automation  Industries  Inc.),  cured  at  79-80°F  for  a minimum  of 
12  hours  under  clamping  pressure.  The  manufacturer  states  a 
20  percent  maximum  usable  strain  at  room  temperature  for  this 
cement. 

The  load  cell  was  a stiff  axial  rod  strain-gauged  type  (Instron 
Model  GR) , with  a ±20,000  pound  full  scale  capacity,  and  ac  sig- 
nal conditioning.  The  output  amplifier  was  capable  of  producing 
full  scale  recorder  pen  motion  for  ranges  as  low  as  500  pounds  at 
rated  accuracy  of  ±-J-  percent. 

The  plate  tests  were  performed  using  the  apparatus  described 
previously.  Alignment  of  the  specimen  and  fixture  was  accomplished 
by  replacing  the  flat-faced  loading  rod  (see  Fig.  5)  with  a sharp- 
pointed  conical  alignment  rod.  The  moving  crossbeam  was  lowered 
until  the  rod  point  indented  the  lower  crossbeam  of  the  machine, 
marking  the  center  of  the  load  axis.  The  support  fixture  was 
aligned  to  this  center  and  held  in  place  by  two  steel  pins.  The 
three  LVDT's  were  then  positioned  and  held  by  double-coated  mask- 
ing tape  to  the  lower  crossbeam. 
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The  LVDT’s  were  calibrated  at  this  point  by  screwing  a thick 
circular  flat  plate  onto  the  end  of  the  alignment  rod.  This  cali- 
bration plate  was  small  enough  in  diameter  to  fit  within  the 
5.22-inch  support  plate  opening.  While  moving  the  crosshead 
downward  by  known  amounts  of  displacement,  the  signal  balancing 
unit  and  strip  chart  recorder  amplifiers  were  adjusted  to  effect 
the  calibration. 

The  plate  specimen  was  then  installed  and  the  strain  gauge 
leads  connected.  Each  plate  had  a small  center  punch  mark  on  its 
top  face,  which  was  used  to  align  the  plate  by  matching  that  mark 
with  the  point  of  the  alignment  rod.  The  plate  was  then  held  in 
position  by  three  sets  of  bolts  and  washers  around  its  edge.  The 
load  cell  circuit  was  zeroed  and  calibrated  (500  pound  internally 
simulated  load),  and  the  strain  gauge  and  LVDT  circuits  were  checked 
and  zeroed.  The  strain  gauge  scanner  unit  was  set  to  scan  through 
all  the  gauges  in  sequence  and  repeat  continuously,  at  the  nominal 
rate  of  0.2  sec  per  gauge.  The  conical  alignment  rod  was  removed 
from  the  load  cell  and  replaced  by  the  flat-faced  loading  rod.  The 
crosshead  was  lowered  to  bring  the  loading  rod  in  contact  with  the 
plate  surface,  a five-pound  preload  was  applied  to  hold  the  plate 
secure,  and  then  the  holddown  bolts  and  washers  were  removed. 

The  crosshead  motion,  the  strain  gauge  scanner/recorder , tbe 
LVDT  strip  chart  recorder,  and  the  load  cell  chart  recorder  were 
started  simultaneously.  From  that  moment  until  the  end  of  the 
test,  the  operation  was  automatic  and  monitored  continuously.  The 
displacements,  strains,  and  load  values  were  all  related  by  the 
time  parameter,  since  the  various  chart  speeds  and  strain  gauge 
scanner  speed  were  all  accurately  measured  in  advance. 
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Material  Properties  Tests 


Eight  tension  coupons  and  four  compression  rods  were  tested 
for  both  the  and  inch  plate  material.  The  tension  coupon 
tests  were  also  performed  in  the  Instron  machine  (Fig.  4),  in  the 
upper  part  above  the  moving  crossbeam.  The  load  cell  was  attached 
through  the  upper  crossbeam,  with  a universal- joint  coupling  into 
an  upper  grip,  which  held  the  upper  end  of  the  tension  coupon. 

These  items  are  not  shown  here.  The  coupon  was  held  from  below  by 
a lower  grip,  shown  in  Fig.  4,  rigidly  attached  to  the  moving  cross- 
beam. The  tension  was  applied  by  the  downward  moving  crossbeam. 

The  load  was  measured  using  the  Instron  load  cell,  while  the  strain 
was  measured  by  both  the  strain  gauges  and  a clip-on  extensometer . 

The  strain  gauges  were  from  the  same  lot  as  those  used  on  the 
plates,  with  gauges  bonded  to  opposite  faces  of  the  tension  coupons. 

A typical  tension  coupon  with  gauges  installed  is  shown  in  Fig.  9. 
These  strain  gauges  were  applied  using  the  same  techniques  and  mate- 
rials described  previously  for  the  gauges  on  the  plate  specimens. 

For  each  plate  thickness,  the  four  coupons  numbered  1,  2,  8,  and  9 
had  the  dual-element  gauges  while  the  four  numbered  3,  6,  11,  and 
12  had  the  single-element  gauges. 

The  extensometer  (Instron  Model  G-51-11)  was  a strain-gauge 
element  type  that  contacted  one  face  of  the  coupons  with  dual  knife 
edges  held  by  spring  clips.  This  unit  had  a one-inch  nominal  gauge 
length,  a maximum  extension  of  0.1  inch,  and  an  accuracy  of 
percent  (manufacturer’s  calibration  data).  The  extensometer 
signal  was  recorded  through  an  amplifier  as  chart  drive  motion  on 
the  Instron  recorder.  The  extensometer  was  used  to  check  out  the 
accuracy  of  the  strain  gauges  in  the  high  strain  ranges  (1  to 
10  percent),  for  all  eight  of  the  -q— inch  thick  coupons  and  for 
two  (one  with  single-element  and  one  with  dual-element  strain  gauges) 
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of  the  £-inch  thick  coupons.  The  net  strain  gauge  error,  after 
data  processing,  is  discussed  in  Appendix  B. 

The  tension  coupons  were  tested  to  about  10  percent  strain 

in  the  Instron  testing  machine  with  a universal-joint  upper  grip 

and  a fixed  lower  grip.  The  crosshead  speed  was  set  at  values 

that  varied  among  the  tests  from  0.002  to  0.020  inch  per  minute. 

For  the  tension  coupons  used  here,  the  nominal  strain  rate  in  the 

gauge  section  was  estimated  to  be  in  the  range  of  r00  to  6000 

-6 

microstrain  (10  inches/inch)  per  minute. 

The  load  and  extension  were  recorded  simultaneously  on  the 
Instron  chart  recorder,  with  the  load  cell  driving  the  pen  motion, 
and  the  extensometer  driving  the  chart  motion,  producing  x-y  type 
records.  The  strain  gauges  on  the  tension  coupons  were  sampled  by 
scanning  through  the  set  of  gauges  at  different  levels  of  indicated 
load.  The  strain  gauge  scanner-recorder  unit  was  set  to  scan  once 
through  the  set  of  2 or  4 gauges,  at  the  rate  of  0.1  sec  per 
gauge,  therefore  taking  0.2  and  0.4  sec  for  each  scan.  The 
corresponding  maximum  changes  in  strain  during  a scan  were  estimated 
to  be  20  and  40  microstrain. 

The  strain  gauge  data  were  processed  by  mathematically  compen- 
sating for  the  known  errors  due  to  transverse  sensitivity  and  bridge 
nonlinearity,  as  described  in  Appendix  B. 

For  the  compression  tests,  the  load  cell  was  attached  to  the 
lower  crossbeam,  with  a hardened  steel  compression  plate  on  the 
load  cell  upper  end.  Another  compression  plate  was  attached  to  the 
underside  of  the  moving  crosshead.  The  compression  coupons  were 
placed  on  the  load  cell  plate  and  were  compressed  by  the  downward 
moving  crosshead  plate  at  a constant  speed  of  0.002  inch  per  minute. 
This  produced  a strain  rate  of  2000  microstrain  per  minute  on  the 
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one-inch  long  rod  specimens.  These  tests  were  stopped  at  about 
2 percent  strain  because  calculations  and  exploratory  tests  in- 
dicated a tendency  for  plastic  buckling  at  higher  strains. 

The  Instron  recorder  was  operated  in  a strip  chart  mode,  with 
the  load  plotting  across  as  the  chart  moved  down  at  a known  constant 
speed.  The  compression  was  recorded  indirectly  by  relating  the 
chart  motion  to  the  crosshead  motion  and  adjusting  for  the  cali- 
brated deformations  of  the  machine. 
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3 . DATA 


Material  Properties 

The  tension  and  compression  coupon  tests  provided  data  to  de- 
termine the  true  stress-strain  curves  for  the  plate  material.  The 
corrected  strain  gauge  data  were  used  for  this  purpose,  and  the 
details  of  the  strain  gauge  data  processing  are  given  in  Appendix  B. 

A typical  load  versus  extension  record  is  shown  in  Fig.  10, 
which  indicates  that  the  material  fad  a nonsmooth  tension  stress- 
strain  curve  at  high  plastic  strains . This  was  most  probably  due 
to  a nonuniform  distribution  of  yielding  known  as  the  Portevin- 
le  Chatelier  effect  (Ref.  20  and,  for  a more  recent  example, 

Ref.  21),  in  which  incremental  plastic  strains  are  confined  to  con- 
centrated yield  bands,  similar  to  Liiders  bands  in  steel.  A smooth 
average  curve  could  be  faired  through  this  trace,  with  little  dif- 
ficulty. The  digitalized  strain  gauge  data,  however,  were  recorded  at 
nonuniform  discrete  intervals,  which  appeared  to  add  a certain  amount 
of  random  scatter,  but  it  was  still  possible  for  an  accurate  curve 
to  be  faired  through  the  digital  data,  as  shown  in  Fig.  11. 

In  this  way,  stress-strain  curves  were  constructed  for  each 
tension  coupon,  and  from  these,  true  stress  versus  average  strain 
curves  were  calculated  by  averaging  the  individual  values  of  cor- 
rected strain  gauge  data  at  intervals  of  true  stress.  The  method 
of  calculating  true  stresses  is  given  in  Appendix  C.  These  curves 
are  shown  in  Figs.  12  and  13  for  the  and  £-inch  plates,  and 
an  analysis  of  the  variations  in  the  coupon  data  are  given  in 
Tables  2 and  3.  The  tables  indicate  that  the  plate  material  was 
essentially  isotropic  and  homogeneous,  since  the  differences  be- 
tween the  tension  properties  in  the  longitudinal  and  transverse 
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TABLE  3 SUMMARY  OF  COUPON  TESTS  - J-INCH  PLATE 
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directions  were  small  and  the  deviations  of  individual  coupons  from 
the  average  were  generally  also  small.  One  exception  was  the  elas- 
tic moduli.  These  measurements  were  made  more  difficult  by  the 
small  elastic  strain  limit  of  the  material,  so  that  the  variability 
was  probably  caused  more  by  instrumentation  and  data  interpretation 
errors  than  by  any  actual  variation  in  the  moduli. 

Note  that  only  2 compression  tests  are  reported,  both  for  the 
inch  plate.  Data  from  the  other  compression  tests  were  judged 
to  be  invalid  because  the  elastic  slopes  of  those  stress-strain 
curves  did  not  agree  sufficiently  with  the  known  elastic  modulus 
for  this  material.  Those  specimens  were  too  stiff  for  satisfactory 
use  in  the  testing  machine.  The  specimen  length  had  to  be  kept 
short  to  prevent  plastic  buckling,  but  this  caused  the  specimen 
stiffness  to  exceed  the  machine  stiffness.  Large  errors  in  the 
calculated  specimen  strains  resulted,  because  of  apparently  uncon- 
trollable variations  in  machine  stiffness  from  test-to-test.  This 
problem  could  have  been  avoided  by  repeating  the  tests  with  a one- 
inch  gauge  length  compres some ter  to  directly  measure  the  deforma- 
tions on  the  specimens.  It  was  decided  to  abandon  further  com- 
pression tests,  however,  since  the  two  acceptable  tests  confirmed 
the  published  information  that  the  compression  and  tension  stress- 
strain  curves  were  nearly  identical,  as  shown  in  Fig.  12. 

A mathematical  representation  of  the  stress-strain  curves  of 
the  type 

e = (a/E)  + aan  (1) 

was  required  for  use  in  the  theoretical  prediction  of  the  plate 
behavior,  and  therefore  such  curves  were  fitted  through  the  test 
data.  The  first  term  on  the  right  hand  side  of  Eq.  (1)  is  the 
elastic  strain  and  the  second  the  plastic  strain.  The  approxi- 
mating curves  plotted  in  Figs.  12  and  13  represented  the  best 
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compromises  to  fit  over  the  whole  range  of  strains.  Apparently, 
the  plastic  strains  did  not  fit  the  second  term  of  Eq.  (1)  very 
well  over  the  range.  A least  squares  curve-fitting  analysis 
showed  that  a slightly  better  fit  would  result  with  a formula 
using  aencr  for  the  plastic  strain  term.  However,  theoretical 
analysis  requires  that  the  stress-strain  law  be  in  the  form  of 
Eq.  (1). 

A better  theoretical  approach  would  be  to  approximate  the 
stress-strain  curve  by  piecewise  linear  segments.  This  would 
allow  the  representation  to  be  as  accurate  as  desired,  by  varying 
the  number  of  points  connected  by  the  straightline  segments. 

Plate  Behavior 

Test  data  from  the  three  plates  are  discussed  in  this  section. 
For  each  plate,  a basic  set  of  data  includes: 

• Data  histories  - deflections  and  strains  plotted 
versus  load,  at  regular  intervals  up  to  the  maxi- 
mum test  load 

• Strain  distributions  - strains  on  lower  (tension) 
surface  plotted  versus  radial  position  at  several 
load  levels 

• Symmetry  check  - radial  strains  from  three  gauges 
equally  spaced  around  a 1.0-inch  radius  circle, 
on  upper  and  lower  surfaces,  plotted  versus  load. 

Also  to  be  discussed  are  the  corresponding  theoretical  pre- 
dictions of  the  Grumman  finite  element  computer  code  AXSHEL  (Axi- 
symmetric  Shell  Analysis)  of  the  PLANS  (Plastic  and  Large  deflec- 
tion Analysis  of  Structures  System).  The  version  of  AXSHEL  used 
here  is  capable  of  treating  the  combined  nonlinearities  associated 
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with  elastic-plastic  material  behavior,  and  large  deformations 
for  bending  and  membrane  behavior.  It  does  not  include  trans- 
verse shear  stiffness  nor  local  (three  dimensional)  stress  dis- 
tributions at  the  loaded  areas.  A very  brief  description  of 
the  general  theoretical  method,  and  the  specific  finite  element 
idealizations  used  here  is  given  in  Appendix  D.  A more  complete 
description  can  be  found  in  Refs.  22,  23,  and  24. 

Data  for  the  |-inch  thick  plate  No.  1A125,  having  a supported 
diameter-to-thickness  ratio  of  D/h  s;  40.6  is  given  first,  in 
Figs.  14  through  19.  This  plate  was  tested  to  a maximum  load  of 
1140  pounds,  developing  a peak  deflection  of  approximately  2.6 
times  the  thickness.  Two  finite  element  theory  models  were  used, 
one  having  14  and  the  other  27  annular  plate  elements  (Fig.  D-l) . 
load  increments  of  2,9  pounds  were  applied  to  both  models. 

The  loading  was  modeled  by  assuming  a uniform  contact  pressure 
between  the  load  rod  and  plate.  It  is  known  that  the  contact  pres- 
sure between  a rigid  rod  and  an  elastic  surface  is  highly  nonuni- 
form with  more  contact  pressure  at  the  edge  than  at  the  center  of 
the  rod.  However,  it  was  decided  to  try  to  use  the  uniform  con- 
tact pressure  as  a simple  first  approximation  to  the  more  complex 
actual  distribution,  since  the  loaded  area  was  small  compared  to 
the  whole  plate;  only  the  region  near  the  load  was  expected  to  be 
affected. 

The  deflection  data  for  this  plate  are  shown  in  Fig.  14.  The 
main  feature  to  note  is  the  distinctly  nonlinear  behavior,  with 
plasticity  at  first  causing  reduced  stiffness,  followed  by  a re- 
stiffening of  the  plate  at  higher  loads  as  large  deflections 
(w  /h  > 1)  are  developed.  The  measured  motion  of  the  load  rod, 
in  contact  with  the  upper  surface,  is  plotted  in  Fig.  14  as  the 
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circular  data  points.  The  square  data  points  show  the  data  from 
the  central  LVDT,  which  was  in  contact  with  the  lower  surface  at 
the  center.  That  LVDT  had  a useful  range  at  about  0.22  inch, 
so  only  the  load  rod  data  are  available  after  that  point. 

The  predicted  displacements  curves  are  drawn  for  the  loca- 
tion R = 0.215  inch,  under  the  edge  of  the  loading  rod,  as  well 
as  for  the  locations  of  the  three  LVDT  transducers.  The  predic- 
tions agree  very  well  with  the  data  for  the  two  outer  positions. 
However,  both  the  load  rod  motion  and  the  central  LVDT  data  cor- 
relate more  closely  with  the  prediction  for  the  motion  of  the  load 
rod  edge  at  R = 0.215  inch,  than  with  the  predicted  central  de- 
flection at  R = 0.  It  could  be  expected  that  the  load  rod  would 
tend  to  contact  the  plate  at  its  edge  circle  as  the  deformation 
progressed  and  the  data  seem  to  confirm  that.  However,  the  mea- 
sured central  deflection  being  nearly  equal  to  the  rod  edge  de- 
flection indj' ites  that  the  plate  remained  nearly  flat  under  the 
load  rod  (at  least  up  to  0.22  inch  deflection). 

This  difference  between  the  measured  and  predicted  central 
deflections  was  probably  caused  by  inaccuracy  of  the  assumed 
contact  pressure  distribution.  Away  from  the  rod  edge,  the  plate 
behavior  cannot  depend  significantly  upon  the  load  distribution, 
but  only  upon  the  total  load  level. 

Figure  14  shows  that  nearly  doubling  the  number  of  elements 
in  the  theoretical  model,  from  14  to  27,  produced  a small  in- 
crease in  predictive  accuracy  for  deflections,  resulting  in  excel- 
lent correlation  with  the  test  data,  but  the  computing  cost  was 
nearly  doubled.  Probably,  less  than  14  elements  would  have  been 
sufficient  to  predict  the  deformations,  if  an  error  of,  say,  five 
percent  would  have  been  tolerable. 
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The  processed  strain  gauge  data,  as  a function  of  applied  load, 
are  shown  in  Figs.  15  and  16.  The  lower  (tension)  surface  data  are 
replotted  in  Figs.  17  and  18,  as  strain  distributions  along  a 
radius,  for  several  representative  load  levels.  The  upper  surface 
circumferential  gauge  and  both  radial  gauges  at  R = 0.51  inch 
failed  after  300  pounds  load,  perhaps  because  they  were  too  close 
to  the  loading  rod.  There  was  no  strain  gauge  at  R = 0 on  the 
top  surface,  but  the  theoretical  curve  is  given  for  that  location 
for  reference  purposes. 

The  most  noticeable  features  of  the  strain  data  are  chat  the 
circumferential  strains  were  greatly  dominant  over  most  of  the 
plate,  while  the  radial  strain  component  drops  rapidly  with  radius, 
being  significantly  lerge  only  in  the  region  under  the  load  circle. 
Apparently,  the  flow  mechanism  is  primarily  circumferential,  but  a 
radial  "plastic  hinge"  effect  is  important  near  the  plate  center. 
This  is  qualitatively  consistent  with  previously  reported  results 
for  mild  steel  simply  supported  circular  plates  (Refs.  5 and  14). 

The  correlation  between  measured  and  predicted  strains  is  good 
only  on  the  outer  half  of  the  plate  away  from  the  load  rod,  where 
the  strains  are  relatively  low.  The  theory  greatly  overpredicted 
the  tension  strain  near  the  plate  center,  with  the  error  increas- 
ing with  greater  load  levels.  At  the  l^ad  level  of  1017  pounds 
the  predicted  strain  at  the  center  was  almost  twice  the  measured 
value  (Fig.  17).  The  predictive  accuracy  improved  with  radius, 
as  the  strain  level  dropped.  In  general,  the  overprediction  in- 
creased with  increasing  strain. 

This  large  error  in  strain  predictions  at  the  plate  center 
was  most  likely  caused  by  the  assumed  distribution  of  contact  pres- 
sure, as  previously  described.  Since  the  strain  is  usually  more 
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sensitive  to  local  conditions,  it  is  to  be  expected  that  the  strain 
predictions  under  the  load  would  be  more  in  error  than  the  deflec- 
tions . 

Another  possible  cause  for  the  strain  prediction  error  might 
be  that  the  theory  does  not  include  the  effects  of  transverse  shear 
or  of  local  contact  stresses.  Since  the  latter  will  be  generally 
compressive  in  the  region  of  the  plate  under  the  loading  rod,  it 
is  to  be  expected  that  the  measured  tensile  strains  will  be  reduced 
at  the  plate  center.  Timoshenko  (Ref.  25)  discusses  this  case  for 
purely  elastic  behavior,  but  provides  a formula  for  the  peak  ten- 
sion strain  only  for  the  limiting  case  of  a uniformly  loaded  circle 
of  very  small  size  compared  to  the  thickness.  Since  the  case  at 
hand  has  a loaded  circle  diameter  2c  = 0.415  inch  and  a thickness 
h = 0.128  inch,  the  ratio  2c/h  = 3.24  is  outside  the  range  of 
applicability  of  Timoshenko’s  formula.  Anderson  and  Shield  (Ref.  26) 
examined  the  problem  using  perfect  plasticity  and  the  Tresca  yield 
condition.  They  found  that  the  local  behavior  under  the  load 
governed  the  plastic  collapse  load  for  2c/h  < 1,  while  plate 
theory  was  sufficiently  accurate  to  predict  the  collapse  load  for 
2c/h  > 3.  The  three  plates  tested  here  have  values  of  1.4  < 

2c/h  <3.2  (listed  in  Table  4),  that  are  in  the  range  where  the 
contact  stresses  can  be  expected  to  have  an  effect  on  the  over-all 
plate  behavior;  for  example,  on  the  deflections. 

The  strain  predictions  using  a 27-element  theoretical  model 
(Appendix  D)  are  plotted  for  two  load  levels,  as  the  dashed  curves, 
Figs.  17  and  18.  As  with  the  deflections,  they  show  very  little 
difference  from  the  14-element  model  and  a still  fewer  number  of 
elements  would  have  been  sufficient. 

At  the  plate  center,  the  radial  and  circumferential  strains 
are  identical,  as  required  by  the  symmetry,  so  that  a single  strain 


Table  4 LOAD  CONCENTRATION  PARAMETER 


Plate  Specimen 


1A125 

4A125 

4A250 

2c 

0.415 

0.375 

0.375 

h 

0.1283 

0.1286 

0.2615 

2c/h 

3.2 

2.9 

1.4 

gauge  was  used  for  both  strain  components  at  that  location.  The 
actual  symmetry  of  the  plate  test  was  checked  by  comparing  the 
strain  data  from  the  radial  gauges  installed  along  a 2-inch  diame- 
ter circle  at  120  degree  intervals.  These  strain  data  are  shown 
in  Fig.  19.  Evidently  the  symmetry  on  the  bottom  surface  shows  some 
variance  among  the  three  gauges.  This  leads  to  a degree  of  uncertainty 
about  the  accuracy  of  the  test  data  for  this  plate,  because  the  varia- 
tion is  on  the  surface  having  the  largest  strains. 

Plate  test  4A125  was  essentially  a duplicate  of  the  first  test, 
except  that  the  load  rod  diameter  was  reduced  to  0.375  inch  to 
increase  the  distance  between  the  rod  edge  and  the  innermost  strain 
gauges,  in  hopes  of  preventing  any  premature  gauge  failures. 

The  test  data  and  predictions  of  the  finite  element  theory 
are  shown  in  Figs.  20  through  24.  Test  data  for  this  plate  were 
qualitatively  similar  to  that  for  the  first  plate,  with  slight 
increases  in  the  central  deflection  and  strains  near  the  load  that 
are  attributed  to  the  greater  load  concentration.  Two  sets  of 
theoretical  predictions  were  made  for  this  plate.  One  assumed  ?. 
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uniformly  distributed  contact  pressure  under  the  loading  rod, 
before,  and  tht  other  assumed  a line  load  concentrated  along  the 
edge  circle  of  the  loading  rod.  The  circular  line  load  case  was 
added  because  ttie  data  from  the  previous  plate  indicated  that  the 
assumption  of  uniform  contact  pressure  was  inadequate.  The  actual 
conLact  pressure  was  expected  to  be  somewhere  between  the  two 
limiting  cases  of  uniform  pressure  and  circular  line  load. 

The  load  versus  deflection  data  in  Fig.  20  show  only  the  mea- 
sured motion  of  the  load  rod,  because  an  electrical  circuit  failure 
between  the  three  LVDT  displacement  transducers  and  recorder  elimi- 
nated their  records. 

Figure  20  shows  various  types  of  theoretical  predictions  for 
deflection,  generated  by  AXSHEL  as  it  was  being  developed  during 
this  test  project.  The  solid  curves  are  for  the  predictions  that 
assumed  a uniformly  distributed  contact  pressure  under  the  loading 
rod.  At  first,  the  theory  contained  only  the  elastic-plastic  mate- 
rial behavior  with  small  deflections  and  strains  (the  lower  curve). 
But  the  test  data  showed  that  after  plasticity  began,  the  plates 
were  stiffer  than  predicted  by  that  restricted  theory.  The  small 
deflection  plastic  bending  theory  overpredicted  the  deformations 
and  strains,  starting  with  a central  deformation  of  about  40  per- 
cent of  the  plate  thickness  (w  /h  « 0.4).  This  indicates  that 
the  membrane  strains  and  stiffness  became  significant  at  that  point. 
The  prediction  seemed  to  be  approaching  a collapse  load  — a zero- 
stiffness condition  defined  as  a large  deflection  increase  for  a 
very  small  load  increase  — at  about  500-600  pounds.  The  test 
data,  however,  showed  no  collapse  trend,  but  an  indication  to  re- 
gain stiffness  continually  after  w^/h  = 0.7.  This  trend  was  still 
progressing  when  the  test  load  was  arbitrarily  stopped  at  about 
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The  linear  elastic  (small  deformation)  case  is  shown  for  ref- 
erence purposes.  Diverging  from  it  is  the  elastic  large  deflection 
prediction  that  incorporates  the  effect  of  the  membrane  stresses, 
generated  from  large  deformations,  on  the  bending  stiffness  of  the 
plate.  Both  of  these  grossly  underpredict  the  deformations,  as 
expected,  because  of  the  neglect  of  plastic  behavior. 

Finally,  the  elastic-plastic  large  deflection  theory  is  plotted 
that  includes  both  bending  and  membrane  stiffnesses.  The  solid  curve 
is  for  the  uniform  contact  pressure  and  the  dashed  curve  is  for  the 
concentrated  line  load  along  the  loading  rod  circle.  Up  to  about 
400  pounds,  the  predictions  for  the  two  loadings  are  the  same,  with 
an  excellent  agreement  with  the  test  data.  At  larger  loads  the  theo- 
retical deflections  separate,  with  the  circular  line  loading  case  pro- 
ducing an  excellent  prediction  of  the  test  data. 

The  strain  data  are  shown  in  Figs.  21,  22,  and  23,  along  with 
the  predictions  for  the  two  theoretical  loading  cases.  As  expected, 
the  two  theoretical  results  are  nearly  the  same  away  from  the  loaded 
area,  but  diverge  greatly  at  the  center,  under  the  load.  The  strain 
data  agree  fairly  well  with  the  theory  for  the  circular  line  load, 
except  on  the  upper  surface  at  R = 0.53  inch,  and  on  the  lower 
surface  at  the  center,  R = 0.  The  measured  strains  at  the  center 
of  the  lower  surface  fall  between  the  two  theoretical  limiting  cases, 
but  closer  to  the  prediction  for  the  circular  line  load.  This  is 
taken  as  a confirmation  that  the  actual  contact  pressure  distribution 
was  somewhat  nonuniform  with  a greater  concentration  of  pressure  at 
the  rod  edge  than  at  the  center.  The  contact  pressure  distribution 
can  be  expected  to  change  with  load,  depending  on  the  local  yielding 
of  material  in  the  contact  region. 

Note  that  the  two  theoretical  loading  cases  produce  signifi- 
cantly divergent  strains  for  R < 0.5  inch.  This  is  equivalent  to 
four  thicknesses  (R/h  < 4)  or  2.7  times  the  radius  of  the  loading  rod. 
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Circumferential  Strain  History,  Plate  4A125 
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RADIUS  R,  INCHES 

Fig.  23.  Radial  Strain  Distributions,  Lower  Surface, 
Plate  iJA125AL 


Of  interest  is  the  symmetry  of  the  strain  gauge  data  in 
Fig.  24.  This  time,  the  strain  gauge  data  indicates  a good  uni- 
formity around  the  circumference  at  R = 1 inch,  and  it  can  be 
concluded  that  the  symmetry  for  this  plate  was  satisfactory. 

The  strain  data  for  this  plate  were  examined  to  determine  the 
relative  magnitude  of  the  membrane  (mid-surface)  strains  and  bend- 
ing strains.  The  motivation  was  to  see  if  the  cause  of  the  pre- 
dictive error  could  be  traced  to  a particular  mechanism,  either 
bending  or  membrane,  so  that  improvements  in  the  theory  could  be 
specifically  directed.  This  was  done  simply  by  assuming  a linear 
variation  of  strain  through  the  thickness  of  the  plate.  This  is 
consistent  with  thin  plate  theory,  in  which  the  thickness  contrac- 
tion due  to  in-plane  strains  is  not  significant.  The  membrane  and 
bending  strains  are  then 
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The  circumferential  membrane  and  bending  strain  distributions, 
calculated  from  the  test  data,  are  plotted  in  Fig.  25  along  with 
the  corresponding  theoretical  predictions  of  the  AXSlffiL  code.  Note 
that  the  predictions  are  generally  more  accurate,  compared  to  the 
test  data,  for  the  membrane  strains  than  for  the  bending  strains. 
This  might  lead  to  the  judgment  that  the  bending  theory  was  pri- 
marily responsible  for  any  lack  of  agreement  between  predicted  and 
measured  strains.  However,  the  bending  strain  levels  were  gene- 
rally greater  than  the  membrane  strain  levels,  and  the  predictive 
error  might  be  caused  by  a general  overprediction  of  plastic 
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strains.  Such  a plastic  strain  error  would  be  greater  for  greater 
strain  level,  regardless  of  the  mechanism,  and  could  have  been  the 
main  cause  of  the  discrepancies . 

If  a comparison  is  made  purely  on  the  basis  of  strain  level, 
it  appears  that  the  predictive  error  is  about  the  same  for  the 
bending  and  membrane  strains.  However,  the  membrane  strains  at 
the  gauge  sites  are  not  large  enough  to  make  this  conclusion  defi- 
nite. Therefore,  a clear  judgment  concerning  the  theoretical 
mechanism  requiring  improvement  cannot  be  made  at  this  time. 

Some  interesting  behavior  can  be  seen  by  comparing  Figs.  25 
and  22,  and  recalling  that  the  lower  surface  strain  at  any  point 
is  the  sum  of  the  membrane  and  bending  strains.  It  is  known 
(Ref.  25,  p.  415,  for  example)  that  the  membrane  strains  in  an 
initially  flat  elastic  plate  are  approximately  proportional  to  the 
square  of  the  central  deflection,  during  the  early  stages  of  the 
large  deflection.  Note  that,  as  expected,  the  membrane  strains 
were  extremely  small  at  low  loads,  but  grew  very  rapidly  with  in- 
creasing load  to  become  significant  over  the  entire  plate.  The 
greatest  membrane  strains  were  in  the  central  loaded  region.  In 
the  outer  half  of  the  plate  the  membrane  strains  are  not  large, 
although  they  are  significant  compared  to  the  bending  strains. 

Near  the  outer  edge,  the  compressive  circumferential  membrane 
strain  was  greater  than  the  tensile  bending  strain  and  caused  a 
strain  reversal  on  the  lower  surface  as  the  load  increased  (Figs.  21 
and  22).  Since  the  circumferential  membrane  strain  was  compressive 
outside  a circle  of  approximate  radius  R = 1 inch,  and  since  the 
membrane  strain  levels  were  increasing  more  rapidly  than  the  bend- 
ing levels,  it  seems  safe  to  conclude  that  the  strain  reversal  on 
the  lower  surface  would  have  continued  to  spread  inward  to  cover 
more  of  the  plate. 


This  behavior  is  accounted  for  by  the  AXSHEL  theory,  by  un- 
loading along  the  elastic  slope  of  the  stress-strain  curve  at  the 
points  in  question.  Since  the  kinematic  hardening  rule  is  used, 
the  elastic  unloading  may  cover  a stress  range  of  up  to  twice  the 
initial  elastic  limit  in  tension.  This  may  result  in  a strain 
reversal,  as  reported  here.  In  fact,  AXSHEL  will  account  for  sub- 
sequent plasticity  during  the  reversal,  if  the  new  elastic  range 
is  exceeded.  This  would  be  done  by  calculating  the  changes  in 
plastic  strain  from  the  original  stress-strain  law,  Eq.  (1),  but 
using  reversed  increments  of  stress  and  strain.  This  feature  is 
discussed  in  more  detail  in  Refs.  22,  23,  and  24. 

The  test  data  and  AXSHEL  predictions  for  the  £”inch  thick 
plate  4A250,  with  D/h  = 20.0,  are  shown  in  Figs.  26  through  31. 
The  test  was  conducted  to  a maximum  load  of  4270  pounds,  with  a 
peak  deflection  approximately  1.6  times  the  thickness.  The  theory 
used  14  annular  elements  and  load  increments  of  5.5  pounds,  with 
both  the  uniform  contact  pressure  and  circular  line  load  cases. 

The  load  versus  deflection  data  are  plotted  in  Fig.  26,  which 
again  shows  that  the  measured  central  deflection  follows  closer  to 
the  predicted  deflections  for  the  concentrated  line  loading  than 
for  the  uniform  contact  pressure.  The  center  LVDT  recorded  defor- 
mations almost  equal  to  those  of  the  load  rod,  up  to  approximately 
wq  = 0.16  inch  (wQ/h  ~ 0*6),  at  which  point  the  central  LVDT 
reached  its  usable  limit.  The  appearance  of  these  curves  has  some 
similarity  to  that  of  the  two  ’—inch  plates,  shown  in  Figs.  14 
and  20,  for  equal  w/h  levels.  In  both  cases  the  initial  soften- 
ing is  followed  by  a gradual  restiffening  caused  by  membrane  action 
under  the  large  deflections.  The  main  difference  is  that  for  the 
thicker  plate  of  Fig.  26,  the  plastic  deformation  begins  at  a 
smaller  value  of  w/h  and  the  membrane  restiffening  occurs  at  a 
larger  value  of  w/h  than  for  the  thinner  plates. 
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The  load  versus  strain  data,  plotted  in  Figs.  27  and  28, 
show  that  the  theory  was  generally  in  good  agreement  with  the 
measured  values,  except  on  the  compression  (top)  surface  near 
the  plate  center.  The  largest  strain  component  was  the  tensile 
circumferential  strain  on  the  lower  surface,  for  which  the  theory 
using  the  circular  line  load  provided  very  good  predictions,  ex- 
cept at  the  center.  This  can  best  be  seen  in  Fig.  29,  where  the 
circumferential  strain  distributions  are  plotted.  The  central 
strain  measurement  is  between  the  predicted  values  for  the  two 
limiting  cases  of  load  distribution,  but  closer  to  the  concen- 
trated load  case,  as  for  the  ^-inch  plate. 

The  radial  strain  distribution  on  the  lower  surface  is  plotted 
in  Fig. 39,  which  exhibits  the  very  local  peak  strains  under  the 
load  concentration,  that  is  typical  of  radial  strains  for  this  type 
of  loading.  An  interesting  feature  is  that  the  radial  strain  data 
remained  nearly  constant  over  most  of  the  lower  surface  from  3200 
to  4100  pounds  load,  but  the  theories  predicted  an  increase. 

Note  that  the  effect  of  the  distribution  of  contact  pressure  was 
evident  over  an  increasingly  larger  region  of  the  plate,  as  the 
load  increased. 

The  symmetry  of  strain  was  again  better  on  the  upper  surface 
than  the  lower,  as  indicated  by  Fig.  31. 

The  membrane  and  bending  strains  were  calculated  from  the  test 
data  according  to  Eqs,  (2)  and  (3),  and  are  plotted  in  Figs.  32  and 
33.  These  are  similar  to  those  for  the  4— inch  thick  plates  dis- 
cussed, except  that  the  load  and  strain  levels  are  much  higher.  The 
sum  of  the  bending  and  membrane  strains  (Figs.  32  and  33)  equal  the 
surface  strain  (Fig.  29).  For  this  plate,  the  theoretical  membrane 
strain  predictions  were  lower  in  tension  and  greater  in  compression 
than  the  test  data,  effectively  showing  a compression  shift  relative 
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Co  the  test  data.  The  bending  strains  show  the  opposite  shift  be- 
tween test  and  theory,  however.  The  net  result  was  that  the  two 
predictive  errors  nearly  canceled  each  other,  producing  the  good 
correlation  between  test  data  and  theory  shown  in  Fig.  29.  This 
was  unique  to  the  lower  surface,  since  on  the  upper  surface  the 
bending  strains  reverse  direction  and  then  the  two  errors  were  addi- 
tive, as  can  be  seen  in  Fig.  27,  which  shows  consistently  poorer 
correlation  on  the  top  surface  than  on  the  bottom. 

At  the  highest  load,  the  effect  of  the  theoretical  contact 
pressure  distributions  upon  the  membrane  and  bending  strains  was 
spread  across  the  entire  plate,  even  to  the  outer  edge.  This  was 
surprising,  since  it  seemed  reasonable  to  assume  that  away  from  the 
load  rod,  only  the  load  level  would  determine  the  strains.  Evi- 
dently this  may  not  be  true  for  cases  of  large  nonlinear,  plastic 
deformations.  However,  this  effect  may  have  been  caused  by  the 
numerical  procedures  or  the  finite  element  approximations,  so  that 
using  a greater  number  of  elements  might  have  eliminated  these  dif- 
ferences between  the  two  results. 


4.  CONCLUSIONS 


Test  data  from  these  three  plates  indicated  that  the  qualita- 
tive aspects  of  the  deflection  and  strain  response  were  generally 
similar  for  the  plates  of  diameter-to-thickness  (D/h)  ratios  of 
20.0  and  40.6.  Initially,  at  lower  loads,  the  plates  developed 
a purely  bending  response,  with  plastic  strains  dominating  in  the 
circumferential  direction.  The  radial  plastic  strains  were  large 
in  the  region  of  the  load  concentration  at  the  plate  centers.  As 
the  load  increased,  the  stiffness  (instantaneous  slope  of  load  ver- 
sus deflection,  AP/Awo)  decreased  due  to  the  spreading  of  the 
circumferential  plastic  bending  strains  over  a larger  region  of  the 
plates.  At  some  point,  the  deflections  became  large  enough  to  in- 
duce membrane  strains,  and  brought  membrane  stiffness  into  effect. 
This  occurred  at  a central  deflection  of  approximately  40  percent 
of  the  thickness  (wQ/h  = 0.4).  As  the  load  increased,  the  mem- 
brane strains  increased  more  rapidly  than  the  bending  strains,  and 
the  plate  stiffness  began  to  increase  again  at  approximately 
wQ/h  = 0.8.  This  trend  of  regaining  stiffness  was  continuing  when 
the  tests  were  ended  arbitrarily  at  about  wQ/h  = 2.6  for  the 
thinner  plates  and  wQ/h  = 1.6  for  the  thicker  plate,  although 
the  original  elastic  stiffness  had  not  been  reached.  Evidently, 
as  the  load  increased,  the  bending  stiffness  of  the  plate  decreased 
while  the  membrane  stiffness  increased.  This  can  be  explained  by 
considering  that,  for  any  given  plate,  the  bending  moments  are  pri- 
marily dependent  on  the  load  level  and  material  rigidity,  while  the 
membrane  forces  depend  primarily  on  the  deflection  and  material 
rigidity.  As  the  deflection  induced  by  the  bending  action  accel- 
erates with  load  in  the  presence  of  material  yielding,  the  membrane 
forces  are  induced  and  rise  rapidly  with  load  after  plastic  deforma- 
tion becomes  widespread.  It  is  this  increasing  membrane  action  that 
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accounts  for  the  apparent  restiffening  of  the  plate  as  the  load 
increases . 

The  reason  that  these  opposing  mechanical  trends  can  operate 
simultaneously  in  the  same  body  is  that  the  bending  strains  vary 
linearly  through  the  thickness  and  are  maximum  near  the  surface, 
while  the  membrane  strains  are  uniform  through  the  thickness. 

Thus,  while  the  surface  yielding  is  spreading  and  is  rapidly  losing 
its  resistance  to  the  bending  moments,  the  mid-plane  region  ex- 
periences only  the  smaller  membrane  strains  with  a more  confined 
yielding  area.  Thus,  the  mid-plane  membrane  stiffness  can  still 
be  increasing  while  the  surface  bending  stiffness  drops.  However, 
some  point  must  eventually  be  reached  at  which  the  mid-plane  yield- 
ing will  be  so  widespread  that  the  over-all  membrane  stiffness  must 
also  begin  to  decline.  Then  the  plate  collapse  will  begin,  with  a 
final  flattening  of  the  slope  of  the  load  versus  deflection  curve, 
assuming  that  local  fracture  has  not  already  occurred.  This  final 
yielding  stage  was  not  reached  in  the  tests  reported  here,  although 
it  was  reported  by  Sherbourne  and  Srivastava  (Ref.  18). 

The  finite  element  theory  used  (the  AXSHEL  module  of  the 
Grumman-deve loped  PLANS  system)  produced  excellent  predictions  for 
the  plate  deflections,  and  generally  good  predictions  for  the  strains. 
This  was  true  in  spite  of  the  presence  of  severe  strain  gradients 
caused  by  localized  loading.  The  predicted  values  of  peak  strain 
at  the  plate  center,  assuming  a concentrated  contact  pressure  on 
the  edge  of  the  loading  rod,  were  about  10  to  15  percent  under 
the  measured  strain  values,  for  the  larger  load  levels.  This  dif- 
ference was  attributed  partly  to  the  exact  distribution  of  contact 
pressure  under  the  loading  rod  being  less  concentrated  than  the 
assumed  theoretical  distribution  and  partly  to  the  neglect  of  local 
contact  stresses  in  the  theory.  Further  improvement  might  be  ef- 
fected by  incorporating  an  element  with  three  dimensional  stress 
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states  into  the  AXSHEL  module  for  use  with  load  concentrations  or 
abrupt  changes  in  load  distributions.  This  would  allow  for  a more 
realistic  estimate  of  the  failure  load,  especially  in  cases  where 
local  fracture  may  precede  ductile  collapse. 

Test  data  further  showed  that  the  need  for  including  geometric 
nonlinearity  effects  (large  deformations  and  membrane  forces)  in 
the  theory  was  even  more  important  than  expected.  This  is  because 
these  effects  became  significant  shortly  after  the  onset  of  plastic 
deforma tions. 
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APPENDIX  A 


STRAIN  GAUGE  CONFIGURATION 

The  factory-installed  copper  terminals  of  the  strain  gauges 
can  be  seen  best  in  Fig.  8 (Gauges  15,  17,  and  21),  and  are  con- 
nected to  the  gauge  foil  by  flexible  horizontal  loops  that  are  sup- 
posed to  mechanically  isolate  the  solder  terminals  from  the  gauge 
to  provide  a longer  fatigue  life.  A series  of  exploratory  tests 
was  performed  to  compare  the  peak  usable  strain  level  of  these 
gauges  having  the  factory-installed  terminals  with  that  of  three 
altered  gauge  configurations,  as  follows: 

• Configuration  1:  Unmodified.  Full  carrier  area 

and  factory-attached  copper  jumpers  and  solder  tabs 

• Configuration  2:  Copper  jumpers  removed,  single 

strand  wire  jumper  used 

• Configuration  3:  Same  as  1,  plus  carrier  and  cement 

trimmed  close  around  copper  jumper 

• Configuration  4:  Same  as  2,  plus  carrier  and  cement 

trimmed  close  to  gauge  element 

Configuration  2 is  a common  hand-worked  gauge  installation  technique, 
which  configuration  1 was  designed  to  replace. 

Three  special  2024-0  aluminum  tension  coupons  were  made,  each 
one  having  bonded  to  it  four  gauges.  These  gauges  were  one  each  of 
the  four  configurations  described  above.  The  primary  purpose  of 
these  exploratory  tests  was  to  determine  if  configuration  1 could 
perform  as  well  as  the  standard  method,  configuration  2.  The  other 
configurations  were  added  to  see  if  any  further  improvement  could  be 
a f fee  ted. 
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The  results  of  these  tests  were  that  the  single-strand  j mper, 
configuration  2,  sustained  the  highest  average  failure  strain, 
about  9.5  percent,  with  configuration  1 close  behind  at  8.2  per- 
cent. The  others  performed  relatively  poorly,  in  the  4 to  5 per 
cent  range.  It  was  decided  that  the  convenience  and  labor-saving 
features  of  the  factory-installed  tabs  and  jumpers  more  than  offset 
the  slightly  lower  failure  strain,  and  all  the  strain  gauges  used 
in  this  project  were  of  the  configuration  1 type. 
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APPENDIX  B 


STRAIN  GAUGE  DATA  PROCESSING 


Tension  Tests 

Strain  gauge  data  processing  was  setup  to  correct  the  gauge 
data  for  known  systematic  errors.  Data  were  corrected  by  means  of 
the  formula 

e/e'  = a£  (B-l) 

in  which  e is  the  corrected  strain,  e'  the  indicated  strain, 
and  a and  £ are  the  correction  factors  for  transverse  sensi- 
tivity and  bridge  nonlinearity. 

The  transverse  sensitivity  correction  was  (Ref.  B.l) 

a = (1  - 0.28 5k) / [ 1 + (e^/e')k]  (B-2) 

where  k is  the  manufacturer’s  stated  transverse  sensitivity  factor 
and  is  the  indicated  transverse  strain.  This  formula  is  de- 

rived from  the  fact  that  the  manufacturer's  gauge  specifications 
are  determined  from  calibration  tests  using  steel  bars  having  a 
Poisson’s  ratio  of  0.285.  For  tension  coupons  with  dual  gauge 
sets,  having  transverse  gauges,  the  value  of  a was  recalculated 
for  each  measurement.  For  tension  coupons  with  single  gauges, 
without  transverse  gauges,  the  substitution 

e^/e'  = -|i  (B-3) 

was  made,  where  |i.  was  the  elastic  Poisson's  ratio.  Note  that  in 
the  special  case  of  large  transverse  strains  when  e'/e^  « -k,  the 
value  of  j blows  up,  and  Eq.  (B-2)  cannot  be  used.  In  that  case, 
7 was  interpolated  linearly  from  the  strain  readings  before  and 
after  the  instability.  This  condition  appeared  occasionally  with 
dual  gauge  sets  in  the  plate  tests. 


The  bridge  nonlinearity  correction  was  (Ref.  B.2) 

(3  * 1 + 0.5  Ge#  (B-4) 

where  G is  the  manufacturer's  stated  gauge  factor.  This  correc- 
tion was  required  to  compensate  for  the  slightly  nonlinear  output 
of  the  single  active  arm  type  of  circuit  used  for  these  tests. 

The  extensometer  strain  was  calculated  by  dividing  the  indi- 
cated extension  by  the  initial  gauge  length.  Figure  B-l  shows  the 
remaining  error  of  the  strain  gauge  data,  relative  to  the  extenso- 
meter strain,  after  correcting  by  means  of  Eq.  (B-l).  Note  that 
the  residual  strain  gauge  error  was  negative.  That  is,  the  cor- 
rected gauge  data  were  smaller  than  the  extensometer  data,  and  that 
the  average  error  magnitude  is  slightly  less  than  the  nominal  strain 
level.  Although  the  scatter  in  the  gauge  error  was  large,  the  over- 
all average  error  shows  a systematic  trend  which  could  have  been 
represented  mathematically  and  included  as  another  factor  in 
Eq.  (B-l).  This  was  not  done  primarily  because  the  causes  of  this 
remaining  error  were  unknown  and  it  was  uncertain  how  this  error 
would  be  affected  by  the  bidirectional  strain  states  in  the  circu- 
lar plate  specimens.  Therefore,  it  was  decided  not  to  make  use  .of 
Fig.  B-l  in  correcting  the  strain  gauge  data  for  either  the  tension 
coupons  or  the  plate  specimens,  so  that  the  true  stress-strain  curves 
developed  from  the  tension  data  were  corrected  in  the  same  way  as 
the  strain  data  reported  for  the  plate  tests. 

Plate  Tests 

For  the  plate  tests,  the  strain  gauge  data  were  corrected  for 
transverse  sensitivity  and  bridge  nonlinearity  errors  as  described 
under  tension  tests.  An  additional  correction  was  made  to  compen- 
sate for  the  fact  that  the  radially  and  circumferentially  oriented 
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RESIDUAL  strain  gage  error,  percent 


gauge  sets  had  a slight  rotational  misalignment  with  respect  . j 
the  true  radial  and  circumferential  directions,  as  shown  in 
Fig.  B-2.  This  was  a consequence  of  the  requirement  that  the  dual 
gauge  sets  be  installed  so  that  both  gauges  in  a set  would  be  at 
the  same  radial  distance  from  the  center. 

If  the  gauge  centers  are  apart  by  a distance  26,  then  the 
misalignment  angle  is 

<t>  = arc  sin(5/R)  (B-4) 

Given  the  two  strain  measurements  at  right  angles  to  each  other, 
and  assuming  the  plate  behavior  to  be  perfectly  symmetric  so  that 
any  two  gauges  on  the  same  circle  can  be  considered  to  be  at  the 
same  point,  Mohr's  circle  equations  can  be  used  to  calculate  the 
strains  in  the  exact  radial  and  circumferential  directions. 

In  Fig.  B-2,  the  "radial"  gauge  on  the  left  is  rotated  an 
angle  -<v  from  the  actual  radial  position,  while  the  gauge  on  the 
right  is  an  angle  +1<  from  the  actual  circumferential  position. 
The  Mohr's  circle  for  this  case  is  shown  in  Fig.  B-3.  Note  that 
the  actual  strains,  and  e0 , are  principal  strains,  and  that 

it  does  not  matter  if  the  gauge  positions  are  reversed  (dashed 
lines  in  Fig.  B-3).  The  indicated  strains,  z'  and  are 

- !(cr  + e.)  + i(er  - O cos  2t  (B-5) 

= ?(er  + e(  ) - £(er  - e.)  cos  24  (B-6) 

or 

= 1(1  + cos  2<l)€„  + £(1  - cos  2 <t ) e . (B-7) 

g ' = £(1  + cos  2<t)e  + 1(1  - cos  2<t)er  (B-8) 
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RADIAL  GAUGE  ON  LEFT 


:* 1 z Circle  for  Cirain  Cauri*  Alignment  Error  Analysis 


2 

Since  1 + cos  2:  = 2 cos  i and  1 
become 

2.  . . 2. 

e = cos  ie  + sm  -e 

r r 

2 2 
d'  = cos  + sin  <te 

Solving  for  the  actual  strains,  get 

.2  .2. 

(cos  - - sm  ;)er  - er 

9 2 , 

(cos“^  - sin  l)e,  = e' 
But  from  Eq . (B-4) , 


- cos 

2:  = 2 sin2 1 , 

thes- 

(B-9) 

r 

(B-10) 

2, 
cos  i 

' • 2, 
- e1  sm  i 

(B-ll) 

2 

cos“'t 

/ . 2 
- e s m 1 
r 

(B-12) 

sin2l  = (5/R)2  (B-13) 

cos2l  = 1 - (6/R)2  (B-14) 

and  Eqs.  (B-ll)  and  (B-12)  become 

U - 2(E/R)2Ur  - [1  - (6/R)2]S;  - (6/R)2S;  (B-15) 

(1  - 2(6/R)2 ]e1  = (1  - (5/R)2  \e'  - (6/R)2eR  (B-16) 


Thus  for  any  gauge,  the  correction  factor  is 


1 - (6/R)2 

(6/R)2 

1 - 2 (5/R)2 

1 - 2 (5/R) 2 

( B- 17) 


where  e'  is  the  indicated  strain  of  the  gauge,  and  e,'  is  the 
indicated  strain  of  the  gauge  transverse  to  it.  This  applies  to 
both  the  radial  and  circumferential  gauges. 
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For  the  experiments  described  here,  6 = 0.070  inch  and  the 

minimum  R was  0.5  inch,  so  that  the  maximum  6/R  was  0.14. 

In  that  case  a good  approximation  to  Eq.  (B-17)  can  be  derived  by 

replacing  the  bracketed  terms  by  infinite  series,  and  dropping 
3 

terms  in  (6/R)  or  higher  powers.  The  result  is 

7=1+  (6/R)2  - (6/R)2e^/e/  (B-18) 

Thus,  the  complete  correction  equation  for  the  plate  tests 

was 

e/e'  = (B-19) 

where  */  and  3 were  defined  previously,  and  y comes  from 
Eq.  (B-18). 
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APPENDIX  C 


TRUE  STRESS  CALCULATIONS 


For  the  stress-strain  curves,  the  true  stress  P/A 
over  current  cross  section  area)  was  calculated  from  the 
ing  stress  P/Aq  (load  over  initial  area)  by 


o = 


„ A 
P o 

A A 
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P/A, 


(1  + eT)‘ 


(load 

engineer- 
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where  e is  the  current  strain  transverse  to  the  loading  axis. 
The  is  related  to  the  longitudinal  strain  e,  by 
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where  fi  is  Poisson's  ratio.  This  leads  to 

P/A 

o 

a J 
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(C-2) 


(C-3) 


which  can  be  approximated  by 

o « (1  + 2|Le)P/AQ  (C-4) 

In  the  absence  of  other  data,  it  is  common  practice  to  set  \i  = 

0.5  in  the  plastic  range,  which  gives  a true  stress  of 


o ~ (1  + e)P/AQ  (C-5) 

In  the  tests  described  here,  the  stress-strain  data  from  the 
’-inch  plate  material  were  processed  by  hand.  To  save  time,  for- 
mula (C-5)  was  used.  For  the  r~inch  plate  material,  however, 
the  data  were  processed  by  a digital  computer  program  especially 
made  for  this  purpose.  The  true  stress  values  were  calculated 
using  Eq.  (C-l)  when  transverse  gauges  were  installed,  and  Eq.  (C-3) 
otherwise.  An  effective  value  of  [l  = 0.33  was  used  with  Eq.  (C-3), 
as  this  was  close  to  the  average  effective  value  of  Poisson's  ratio 
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in  the  plastic  range,  calculated  by  the  computer  from  the  trans- 
verse gauge  data. 
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APPENDIX  D 


THEORETICAL  MODEL 

A system  of  computer  codes  is  being  developed  at  the  Grumman 
Aerospace  Corporation,  called  PLANS  (Plastic  Large  deflection  ANaly- 
sis  of  Structures)  which  predicts  the  mechanical  response  of  loaded 
structures  into  the  inelastic  range  of  material  behavior,  using  the 
finite  element  method  of  numerical  structural  analysis.  The  PLANS 
system  treats  elastic-plastic  material  behavior  by  using  the  Von 
Mises  yield  criterion  in  conjunction  with  the  kinematic  hardening 
rule  of  plastic  flow.  Geometric  nonlinearities,  that  is,  altera- 
tions in  structural  characteristics  caused  by  shape  changes  under 
load,  are  also  accounted  for  in  certain  modules  of  PLANS  by  up- 
dating geometric  and  stiffness  properties  in  a piecewise  linear 
fashion.  The  structures  are  modeled  by  the  finite  element  method, 
which  replaces  the  actual  body  by  a number  of  smaller  simply  shaped 
bodies,  or  elements,  connected  at  their  mutually  adjacent  surfaces. 
Since  the  elements  are  relatively  simple  in  shape,  the  mathematical 
model  for  the  deformations  and  strains  in  each  element  type  can  be 
found,  or  approximated,  and  a complex  structure  can  be  approximated 
by  an  assemblage  of  many  simple  elements.  Since  there  is  a separate 
set  of  equations  for  each  element,  the  total  number  of  simultaneous 
equations  to  be  solved  is  generally  large,  requiring  use  of  matrix 
analysis  techniques  in  conjunction  with  a high  speed  digital  computer. 

The  nonlinear  behavior,  caused  by  both  material  plasticity  and 
shape  changes,  is  calculated  by  a series  of  incremental  linear  ap- 
proximations. As  the  load  is  increased  by  a small  increment,  the 
corresponding  increments  of  deformation,"  strain,  stress,  etc.,  are 
calculated  using  the  currently  effective  stiffness  characteristics 
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of  each  element,  as  known  at  the  start  of  each  increment.  Since 
these  effective  element  stiffnesses  are  lagging  slightly  behind  the 
true  stiffnesses,  the  prediction  for  each  increment  will  be  slightly 
in  error.  As  the  previous  incremental  solutions  are  added  up  to 
produce  the  current  total  predicted  state  of  the  structure,  this 
small  incremental  error  will  accumulate  and  increase  as  the  load 
level  increases.  However,  the  accumulated  error  can  be  reduced  by 
using  smaller  load  increments  or  by  performing  certain  corrections. 

Since  the  finite  element  model  assembly  is  an  approximation  to 
the  actual  structure,  the  accuracy  of  the  model  can  also  be  in- 
creased, up  to  a point,  by  using  more  elements  of  smaller  size. 

A description  of  the  PLANS  system,  its  theoretical  basis,  and 
its  operation  is  described  in  Refs.  22,  23,  and  24. 

The  plates  tested  here  were  analyzed  using  the  AXSHEL  module 
of  PLANS.  This  module  handles  axisymmetric  thin  shells  or  thin 
plates,  using  annular  elements.  The  element  behavioral  equations 
use  a shell  theory  which  includes  the  effects  of  bending  moments, 
membrane  forces,  and  large  deflections.  Not  included  in  this  theory 
are  the  effects  of  transverse  shear,  transverse  normal  stresses,  and 
the  local  contact  stresses  in  the  region  of  a load  concentration. 

The  particular  finite  element  models  used  for  each  plate  are 
shown  in  Fig.  D.l.  Either  a 14-annular  element  or  a 27-annular 
element  configuration  was  used,  although  there  were  slight  varia- 
tions among  the  same  basic  configurations  used  for  the  different 
plates.  These  models  differed  from  the  actual  case  primarily  in 
that  the  loading  was  assumed  to  be  either  a uniform  pressure  in- 
side the  loading  circle,  or  a concentrated  load  on  the  circumference 
of  the  circle,  while  the  actual  loading  was  provided  by  a flat-ended 
hard  steel  circular  rod  with  an  unknown  distribution  of  contact 
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pressure.  The  distribution  of -pressure  inside  the  load  circle  was 
expected  to  vary  as  the  plate  deformed  under  load.  Another  devia- 
tion of  the  model  was  that  the  actual  plate  was  slightly  larger  in 
diameter  than  the  support  circle  by  about  0.063  inch,  but  this 
overhang  was  not  accounted  for  in  the  model. 


